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The evolution of the luminosity-temperature-mass 
relations of hot gas in Chandra clusters at 0.4 < z < 1.4. 
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ABSTRACT 

We analyzed the luminosity-temperature-mass of gas {Lx — T ~ Mg) relation for 
sample of galaxy clusters that have been observed by the Chandra satellite. We used 
21 high-redshift clusters (0.4 < z < 1.4). 

We assumed a power-law relation between the X-ray luminosity of galaxy clusters 
and its temperature and redshift Lx ~ (1 -f z)^'^xt T^^-x^ We obtained that for 
an = 0.27 and A = 0.73 universe, Al^t = 1-50 ± 0.23, (3l^t = 2.55 ± 0.07 (for 
68 % confidence level). 

Then, we found the evolution of Mg — T relation is small. We assumed a power- 
law relation in the form Mg ^ (1 + z)"***9^ j^l^MgT a,lso, and we obtained Am^t — 
—0.58 ± 0.13 and PhigT = 1-77 ± 0.16. We also obtained the evolution in Mg — Lx 
relation, we can conclude that such relation has strong evolution for our cosmological 

-/MgLx equation for assuming this 
= 0.73 ±0.15 for a„ = 0.27 



(1 + Z)^"9^ 

i -1.86 ±0.34 and Pm^l, 



parameters. We used Mg 
relation and we found ^m^l 
and A — 0.73 universe. 

In overal, the clusters on big redshifts have much stronger evolution between 
correlations of luminosity, temperature and mass, then such correlations for clusters at 
small redshifts. We can conclude that such strong evolution in Lx—T—Mg correlations 
indicate that in the past the clusters have bigger temperature and higher luminosity. 
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1 INTRODUCTION 

The correlations between different types of X-ray proper- 
ties of galaxy clusters are useful method to obtain informa- 
tion about global properties of clusters. The observation of 
the diffuse. X-ray emitting medium (ICM) of gala:xy clusters 
provides quantities like its mass {Mg), temperature (T), and 
X-ray luminosity {Lx)- The analysis of the scaling relation 
between these physical properties gives considerable insight 
into the physical processes in the ICM. There are difficulties 
for theoretica l predictions o f the evolution of these scaling 
relations (e.g. lNormanll2010l V 

First of all, if concerns with some previous researches 
correlations between X-ray luminos i ty and temperature 
llDavid et al. 1 Il993l . iMarkevitchI Il998l . iMushotzkv fc Scharj 
Il997^ . It is known that there is a correlation between the 
mass of hot intragalactic gas, temperature and X-ray lumi- 
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nosity in clusters. Such a correlation have the small scatter, 
that indicates on the similar formation history for all clus- 
ters. 

However, more detail researches show that the obser- 
vational correlation value is at odds with theoretical pre- 
dictions. The most famous exam ple is the slope Lx ~ T^'^ 
for hot clu sters |Markevitchll 19981 . while the theory provides 
Lx ~ iDavid et al.lll993l " Such differences can be ex- 
plained that heated intergalactic medium play an important 
role non gravitati onal processes a s early heating of massive 
supernovae iCavaliere et al.l (| 19971 ) or radiation cool i ng an d 
the associated active star formation IVoit fc BrvanI (|200lD . 
The evolution of the correlation ratios cause ought to be 
important information that helps to choice between similar 
models. In addition, the correlation value at high redshifts 
is useful for cosmological studies based on the Lx —T — Mg 
evolution of clusters, as they can connect the observed X- 
ray luminosity in ter ms of characteristics of clusters, such as 
temperature or mass iBorgani et"al]|200ll . 

We notice some works on the Lx — T ratio for 
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galaxy clusters at the large red shifts. By Mush otzky with 
so-au t hors in series of paper s (IMushotzkv fc Loewensteini 
19971. iMushotzkv et alj Il996l . IMushotzkv fc Scharj 119971 



Mushotzky et al. 20001 ) analyzed this ratio for large sam- 



ple of distant clusters observed "telescope" ASKA and not 
found t he evolution in Lx — T ratio. How ever, later obser- 
vations (|Arnaud et al.ll2002l . [Novicki et al.i,2002 ') were exis- 
ten ce of such an evolutio n as Lx {z) ~ (l-l-z)^, then A = 1.3- 
l.S. lBorgani et al] (|200lh measurign parameters of 7 clusters 
at 2 > 0.5, concluded about some possible w eak evolution, 
A < 1. A similar conclusion, that made by iHolden et al.l 
l|2002D . measurements for 12 clusters of galaxies. Thus, the 
situation remained be unclear. 

Value Lx — T for clusters has considerable static vari- 
ation. Dispersion in Lx — T for close clusters can be signifi- 
cantly reduced, if in determining the luminosity and temper- 
ature of clusters to exclude the central region ( MarkevitchI 
1 19981 . (|Finoguenov et all l200ll . iFinoguenov et al.ll2007l ')" 



We notice that results from different observational 
methods of mass measurements are not consis tent with one 
another and with the simulat ion results fe.g.. iHorner et ahl 
Il999l . IFinoguenov et al.ll200ll . hereafter FRB). In general, 
X-ray mass estimates are about 80 % lower than the predic- 
tions of hydro-simulations. 

Another recent observational finding is the possible ex- 
istence of a break in the T — Mg relation. By use of resolved 
temperature profile of X-ray clusters observed by ASCA, 
FRB have investigated T — Mg relation in the low-mass end 
and find that M oc r~^, compared to M oc r~^/^ at the 
high mass end. 

In our sample is about 20 clusters a,t z > 0.4 with a 
large exposure time, which allows to determine the temper- 
ature of clusters with an accuracy better 10-15 %. Using 
these Chandra data we were able to accurately determine 
the evolution of the correlation relations between luminos- 
ity, temperature and mass of intracluster gas accumulations 
a.t z > 0.4. 

We derived constraints on the mass-temperature rela- 
tion of gas in galaxy clusters from their observed luminosity- 
temperature relation. We parameterized the Mg—T relation 
as a single power-law, and derive the best-fit values of their 
amplitude and slope from the observed Lx — T relation. 
Also, we derived Lx — Mg relation for our sample, we used 
power-low to describe this relation. 

In Section 2 we show the main steps of the data analysis 
and describe the sample. 

In Section 3 we give the methods to obtain the total 
mass and mass of hot gas. In Section 4, 5 and 6 we show 
the results of our measurements, discussion, and conclusion 
of our work, respectively. 

In this paper we used the cosmological parameters: 
Ho = 73 km s~^Mpc~^ and matter density flrn ~ 0.27. 
Errors are given at the 68% (la) confidence level. 



2 SAMPLE AND DATA PROCESSING 



scribed bv lDickev fc LockmanI (|l99 0h. and coordinates (RA, 
DEC) (which were taken from NEEij). The Chandra obser- 
vations were obtained using the Advanced CCD Imaging 
Spectrometer (ACIS). The list of events (level 1) was pro- 
duced by the Chandra pipeline processing, and then was re- 
processed using the CIAO software package. VFAINT mode 
was used to improve the rejection of cosmic ray events. 

We used the Chandra CIAO software package, version 
4.2, for data reduction, including the latest maps and cal- 
ibration products CalDB 4.5.1 . The main data processing 
steps we re provided using t echniq ue s discussed b y Babvkl 
(l2012ha. iBabyk fc Vayiloval (l2012f) b iBabvk et all (|2012l ')c. 
iBabvk et all (|2012l ')d. iBabvk et all (|2012h e in our previous 
research. 

In brief, we extracted annular spectra centered on the 
X-ray peak with different widths to ensure similar statis- 
tics in the background-subtracted spectra. The X-ray peaks 
were determined with the X-ray images from which all point 
sources were removed. The average of the outermost radius 
is about 1 Mpc. We have generated ARE and RMF files 
using the mkarf and mkrmf command in CIAO tool. Ap- 
proximation of the spectra was carried out in the range 0.5 
- 7.0 keV. 

Temperatures clusters determined by the spectrum col- 
lected within a radius of -R200, this area always contains at 
least 75% of full radiation accumulation. In the approxima- 
tion of the spectra was believed that the interstellar absorp- 
tion in the galaxy corresponds to the thickness of neutral 
hydrogen derived from radio observations. The metallicity of 
intergalactic plasma was considered equal to 0.3 solar value, 
unless the quality of data did not allow to determine it di- 
rectly from the observed spectrum. 

Normalization of the theoretical spectrum, which pro- 
duces the best approximation, can give us the conversion 
factor between the number of photons from the source 
and X-ray luminosity for energy range and extent of is- 
sue EAI = J UeUpdV. This ratio allows to determine the 
emission of intergalactic gas. All luminosity is measured 
within the radius -R200. The central region with a strong 
peak brightness in the center, probably infiuenced by ra- 
diative cooling. In such cases, the central regions (100 kpc) 
were excluded, as the measurement of luminance and spec- 
tral analysis (Markcvitch ct al. 2000). 

The spec tra were analy sed in Xspec (version 12.6) 
environment (I ArnaudI 119961 ). using the MEKAl|3 code 



jKaastra fc Mewe 19931 : incorporating the Fe-L calculations 
of (|Liedahl et al.lll995l l with WAB^ parameter. 

The annular spectra were modeled in order to deter- 
mine the deprojected X-ray temperature and other parame- 
ters, under the assumption of sphe rical symmetry. We have 
used routine metho d DSDEPROJ (|Sanders fc Fabian|[2007l . 
feussell et al.l l200^') to obtain "deprojected" spectra. After 
deprojection, we determined the temperature of clusters and 
build the surface brightness profile for each cluster. 

For each annulus we determined the temperature, kT, 
and parameter norm ~ J n^nudV which is proportional to 



The main observational parameters of the sample galaxy 
clusters are presented in Tab.[T] The table (colums 1-8) con- 
sists of the name of clusters, the redshift, the number of ob- 
servation, the exposition time (before cleaning), the name of 
detectors, and the values of column density, which were de- 



^ http://nedwww.ipac.caltech.edu 

^ MEKAL is the model which describes an emission from hot 
diffuse plasma (ICM) 

Th e WABS is a parameter w hich describe the galactic absorp- 
tion l lDickev fc Lockmanlfiggd) . 
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Table 1. Summary of the Chandra observations (colums 1-8) and result of galaxy clusters (colums 9-11). Colums hst the 
target name, redshift, number of ID observation, net exposure time, detector used, column density and coordinates (taken 
from NED) , also the list of result of galaxy clusters: temperature, luminosities and masses, respectively. 
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the electron, and hydrogen, uh, concentrations. Note, 
that the other parameters of model {Z, z, etc) were fixed. 
Using these parameters, we have built the surface bright- 
ness profiles which fitted by /3 model. The surface brightness 
profiles are measured in the 0.5-7 keV energy band, which 
provides an optimal ratio of the cluster and background fiux 
in the Chandra data. 



3 METHODS 

3.1 Galaxy cluster mass profile 

The gas in a cluster is trapped by the cluster's gravitational 
potential well. If the hot gas is supported by its own pressure 
against gravitational infall, it must obey the equation of hy- 
drostatic equilibrium. The intra-cluster gas in most clusters 
appears to be in approximate hydrostatic equilibrium. 



dP ^ GM{r)p,as{r) 
dr 



(1) 



where P is the pressure of the gas, G is the gravitation con- 
stant, p is the density of the gas and M is the total mass 
inside the sphere of radius r. We used the assumptions of 
spherical symmetry of clusters. Since the density and pres- 
sure of hot gas have very low values, we can use the ideal 
gas law. 



fiirip ' 



(2) 



where n is the mean molecular weight of the gas, and rUp the 
proton mass. Conbining equation ^ and equation ([2]), the 
total mass of cluster as a function of the projected radius is. 



GprUp \ dlnr dlnr 



(3) 



Assuming an isothermal, spherical gas cloud, with tem- 
perature T, equation becomes, 



Ml 



kTr f dlnpg 



Gpnip \ dlnr 



(4) 



where r is the outer radius of the cluster and ^ is constant 
with radius, since the chemical composition of the gas is 
expected to be uniform throughout the cluster. 



3.1.1 Modeling of hot gas emission 

Assuming an isothermal spherical gas cloud in hydrostatic 
equilibrium and assuming th e volume de nsity of galaxies fol- 
lows a so called King profile l|Kindll973 ). the density profile 
of X-ray emitting gas can be approximated by an isothermal 
P model (jCavaliere fc Fusco-Femiaiioll 19761 ). 



Pair) = po 1 + 



-3/2/3 



(5) 



where pg{r) is the density of the gas as a function of the 
cluster's projected radius r, is the core radius and po is 
the density at the cluster's center. The values of /3 and rc are 
obtained by analysis of the X-ray surface brightness profile 
for each cluster. 



3.1.2 X-ray surface brightness profile 

X-ray image analysis allows the determination of the X-ray 
surface brightness profile. Assumption that the hot gas is 
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isothermal leads to a /? model for the cluster's brightness 
profile, 



-3/3 + 1/2 



+ c, 



(6) 



where S{r) is the X-ray brightness as a function of the pro- 
jected cluster's radius r. So, , P and C, are free parameters 
in fitting the model to the X-ray brightness profile. In gen- 
eral the X-ray surface brightness profile is well approximated 
by the /? model. 

There are cases when the standard /3 model leads to 
underestimating the brightness in the central region of the 
clusters. The central X-ray excess is one of the first pieces of 
evide nce of cooling flows in galaxy clusters (jjones fc FormanI 
Il984h. For simpli city, in our study we use a single P model 
( Roncarelli et al. I [2006). In this paper, the X-ray surface 
brightness profile of each cluster in our sample and the fitted 
P model were obtained from the original X-ray images. 



3.2 The gas mass 

The mass of the hot gas in the ICM is determined by the 
integration of the gas density profile (Eq.[5]), over the volume 
of the cluster within a defined radius -R200 , 



A4 = 



pdV. 



Equation ((Tjl can be written as. 



-3/2/3 



dr. 



(7) 



(8) 



The value of po is derived from combining the results 
obtained from fitting a model to the X-ray spectrum, and 
from the best fit parameters obtained by fitting a /? model to 
the X-ray surface brightness profile. The result obtained by 
integrating the density profile over a defined volume of the 
cluster depends on the cosmology assumed because the pro- 



jected radius of the cluster is calculated using the c luster's 
angular distance dA (|Ferramacho fc Blanchardll2007l ). 

The results of calculations are presented in Tab.[T] The 
main source of error in the determination of the gas mass 
comes from the normalization parameter in the spectral fit. 
This parameter is used to compute the gas density. In this 
way, the errors obtained in the gas mass come from fitting 
a single temperature model to the spectrum of a thermally 
inhomogeneous gas. 



4 RESULTS 

Value Lx — T for distant clusters of galaxies, obtained for 
the cosmological parameters Q.a ~ 0.73, ilm ~ 0.27, Ho = 73 
km s^^ Mpc~^, is shown in Fig. [T] Variation in the ratio is 
very small. It is almost entirely due to measurement errors. 
The slope is consistent with the ratio for clusters at low 
redshift. We used normalization parameter on 2; as L ~ (1 -|- 
z)'^lxtt'^lxt and we found that PlxT = 2.55 ± 0.07, and 
AlxT = 1.50 ± 0.23 at confidence level 68 %. 

Now we proceed to the consideration of dependencies, 
which concerns with the mass of hot intergalactic gas. Us- 
ing spherically symmetric clusters, the value of gas mass is 
derived from surface brightness profile. It was decided to 
measure the mass of gas within radius -R200- 

The obtained dependence of Mg, 200 — T are given on 
Fig. [2] We used the normalization parameter as Mg ~ (1 -|- 
z)AM,Trpi3M,T^ as a result we found /Sm^t = 1.77 ± 0.16, 
and ^AfgT = -0.58 ±0.13. 

Relationship between mass and luminosity of gas is 
shown in Fig. [S] This value is useful in assessing the func- 
tions of the distant clusters by their X-ray luminosity. The 
dependence of Mg — Lx has a strong evolutionary relation- 
ship at any values of parameters. We used the normaliza- 
tion parameter as Adg ~ (1 + z)^"'^?^^ L^'^-'g^x ^ so we found 
Pm.Lx = 0.73 ± 0.15, and obtained Am.Lx ~ -1.86 ±0.34 
for A = 0.73, Vim = 0.27 Universe. 

Thus the relationship between X-ray temperature, lu- 
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Deep Cluster Survey in work jRosati et aL 19981) and 




Figure 2. The relation between temperature and mass of gas in 
radius, which correspond 5 =200. The solid line represents Mg—T 
relation for near clusters. 



minosity and mass of gas accumulations at z > 0.4 has a 
strong evolutionary relationship with respect to the correla- 
tions observed at low z. The fact that, as already mentioned, 
the share of baryons in the mass accumulation of full cause 
ought to be universal value, close to the average for the en- 
tire universe (, White et al.,. 1993). but in this case the total 
meiss of baryons can be found as M = Mi,/ fi^. Accordingly, 
the evolution in the ratio Mg — Lx and Mg — T for baryon 
matter and non-baryon cause ought to be the same. 

The observed strong evolution based on Mg — Lx is 
weak, and negative evolution of the luminosity function 
means that the mass function evolves strongly: the spa- 
tial density of clusters of fixed mass at high z is much 
smal ler than in the curent epoch of cosmological parame- 
ters l|Ettori fc Fabianlll999l . lEttori fc FabianlbOOd ). 

Supervisory evolutionary dependence in the ratio Mg — 
T — Lx indicates that the accumulation in the past had a 
great density - they were hotter and had a higher luminosity 
at fixed mass - t hat actually was predict ed scaling theory of 
their formation l|Brvan fc Normanlll998l ). 

For our values of cosmological parameters Q. and A the 
Mg — T relation varies considerably slower than expected 
from theory. Likely explanation for this discrepancy is the 
same as when considerin g the slope of the correlation Lx—T 
l|Frederiksen et al.ll2009l ). 



5 DISCUSSION 

Using the new Chandra cluster temperatures, luminosities 
and mass of X-ray hot intracluster gas, we obtain the Lx — 
T — Mx relations in the 2-11 keV temperature interval 
with greatly reduced scatter. The large redshift sample has 
provided a evidence that clusters do not evolve self-similarly 
in the last ~ 11 Gyr (approximately three quarter of the age 
of Universe) . 

For all data we calculated the PlxT = 2.55 ± 0.07. 
Using the X-ray luminosity function from the ROSAT 



the B rightest Cluster Su rvey in paper ( Ebeling et al.l 
1997) .lBorgani et al.l |l999| ') found Pl^t between 3 and 4. 
iMushotzkv fc Scharj (Il997l'l fouii d Al^t = for all sam- 
ple. In paper of Donahue et al.l l| 19991 ). authors found a 
slightly negative value of A (for fim =0.3, A =0.7, AlxT ~ 
-0.8t?? and = 1.0, Alt = -lA^H, Al^t = 1.5). 
iReichart et al] l| 19991 ) found the Lx — T relation for func- 
tion o f cooling flow- c orrect ed luminosities and tempera- 
tures. iReichart et all (|l999l ) like Pl^t =2.80 ±0.15 and 
AlxT = 0.3 5l?-^t These result s in a good agreement with 
our results. iHorner et al.l l| 19991 ) determined for PlxT ~= 
2.98 ±0.14 and AlvT = 0.02 ± 0.16 (Q.^ = 1.0). 

^ Sadat et al.l l|l998h found a positive ev olution Al^t = 
0.5 ± 0.3 for Qm = 1.0, A = 0.0 and No vicki et al.1 (|2002| ) 
AlxT = 1.1 ± 1.1 for n„, = 1.0, A = 0.0, a nd Al^t = 2.1 
zfc 1.0 for A = 0.7, Qm ~ 0.3, respectively. IVikhlinin et al.l 
(|2002[ ) found AlxT = 0.6±0.3 for A = 0.0, = 1.0 and 
AlxT ~ 1.5 ± 0.3 for a A = 0.7, f2m = 0.3 cosmology. 

Our normalization of the Mg — T relation is also in 
a good agreement with those derive d from recent high - 
resolution numerical simulations (e.g.. iBorgani et al.ll200j ) 
that attempted to model non-gravitational processes in the 
ICM (radiative cooling, star formation, and feedback from 
supernovas) . The primary goal of observationally calibrating 
the Mg — T relation is for use in fitting cosmological models 
to the cluster temperature function. We have presented the 
constraints on the mass-temperature relation of hot gas of 
galaxy clusters. 
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